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1. Introduction 


Gamma-ray bursts (GRBs) are the most luminous objects known; they are luminous enough 
that their positions can potentially be used to help map out large-scale universal structures. GRBs 
are tracers of the stellar matter from which they formed, and they significantly outshine their host 
galaxies so that they can be mapped even when the distribution of underlying galaxies cannot. The 


drawback is that the detection rate is small (95 per year by Swift; []14|]), and the rate at which 
redshifts are measured for the detected bursts is even smaller (roughly 35 per year). 

The angular distribution of GRBs has been studied in detail over the past two decades ([||, [Tj, 


12, 11, 151). Initially, the distribution’s angular isotropy was examined in response to the hypothesis 
that GRBs had Galactic origins. After the cosmological nature of GRBs was established, however, 
the focus of isotropy studies shifted to subsamples having potentially different angular distributions 
(|[^, [T^, 1^, 1^). These studies, originally dependent on GRBs for which redshifts had not 

been measured, became more reliable as the number of GRBs with known redshifts increased, 
but remain limited by small sample sizes. Swift’s compilation of a large number of GRBs having 
known redshifts has reinvigorated large-scale isotropy studies, although the small numbers of GRBs 
in specific redshiff ranges (corresponding fo radial shells) limifs fhese studies to the detection of 
large, pronounced anisotropies. 

We have recently identified a surprisingly large anisofropy suggesfive of clusfering in fhe GRB 
angular disfribufion af around redshiff ([||, ^) in fhe general directions of fhe consfellafions 
of Hercules and Corona Borealis. The scale on which fhe clusfering occurs is disfurbingly large: 
fhe underlying disfribufion of matter suggesfed by fhis clusfer is big enough fo quesfion sfandard 
assumpfions abouf Universal homogeneify and isofropy. Forfunafely, Swiff’s confinued defecfion 
of GRBs makes fhe hypofhesis festable: if fhe anisofropy is affribufable fo sfafisfical sampling, fhen 
fhe clusfer should become less pronounced as more GRBs are defected. 

As of November 2013, the redshifts of 361 GRBs have been measured^ with the sky distribu¬ 
tion shown in Fig.l; most of these GRBs were detected by Swift. This sample represents a 28% 
increase over that used in our previous analysis (283 bursts observed until July 2012). The number 
of GRBs in the 1.6 < z < 2.1 redshift range (||^) (where the cluster resides) has increased from 31 
bursts to 44 bursts, a 42% sample size increase that is large enough to warrant an updated analysis. 


2. Data Analysis 

We analyze the angular distributions in each or our predefined redshiff z bins using fhe k*^ 
nearesf neighbor and fhe boofsfrap poinf radius mefhods. These sfafisfical fesfs are chosen so fhaf 
we can direcfly compare our new resulfs wifh our previous ones. 

The redshiff (radial) bins used in fhis sfudy are chosen fo be fhe same as fhose defined in 
our previous work. In fheory fhe small uncerfainfies associated wifh redshiff measuremenfs allow 
fhe GRB sample fo be easily subdivided info as many redshiff bins as desired. In pracfice, small 
number sfafisfics limif fhe confidence wifh which anisofropies can be defected in each bin. 

We have subdivided fhe fofal sample of 361 bursfs info differenl numbers of redshiff bins rang¬ 
ing from fwo bins (note fhaf a choice of one bin corresponds fo fhe bulk angular GRB disfribufion 

'http://lyra.berkeley.edu/grbox/grbox.php 
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Figure 1: The sky distribution of GRBs with measured redshift. Although the distribution of all GRBs is 
fairly isotropic, extinction causes this sample to miss GRBs near the Galactic plane. 


Table 1: An example of the nearest-neighbor test for four radial groups, with redshift boundaries 
defined in the text. Tabulated numbers represent the KS-test significance that two groups have different 
31®* nearest-neighbor distributions. Boldface type indicates that significant (more than 3(7) differences exist 
between group 2 (1.61 < z < 2.68) and other radial groups. 



Zmin 

gr2 

gr3 

gr4 

grl 

2.68 

0.9999999 

0.942 

0.672 

gr2 

1.61 


0.99904 

0.9999988 

gr3 

0.85 



0.960 


in the plane of the sky) to nine bins. These choices not only allow us to explore the angular charac¬ 
teristics of a variety of radial bins corresponding to redshift intervals, but also allow us to identify 
the redshift range within which any anisotropies lie, should we discover them. The choice of nine 
radial bins provides us with narrow z-bins having the smallest number of bursts per bin (« 40) for 
which we feel we can make reasonable, quantifiable estimates on bulk anisotropies. When choosing 
between 2 and 9 radial divisions, we select the bin sizes that allow us to maintain similar numbers 
of bursts in each radial bin. When we choose a number of bins that does not allow an equal number 
of bursts to be placed in each bin, then we redefine the bin boundaries so that the excess GRBs are 
those with the smallest redshifts; these are subsequently excluded from the analysis. 

The first statistical test we apply to the radially binned distributions is the k**' nearest-neighbor 
statistic; this test looks at the angular separation between each burst and the k**' closest burst to it. 
When this test considers only the nearest neighbor {k = 1), it is sensitive to anisotropies on small 
angular scales corresponding to paired bursts. When looking at widely-separated bursts with large 
k values, the test is sensitive to anisotropies on much larger angular scales. 

We apply the k*^ nearest-neighbor statistic to all burst pairs in each of our radial bins, for 
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the bins in each sample of two to nine radial bins. Most of the times the results are statistically 
consistent with isotropy, regardless of the value of k, the number of radial bins chosen, or the 
particular radial bin being examined. The exceptions arise in the radial bin containing GRBs with 
redshifts of z 2, where anisotropies are observed for medium-sized k values consistent with 
clustering. These anisotropies are most pronounced in the 1.6 < z < 2.1 radial bin, to which this 
clustering appears confined. This result is consistent with our previous findings for a smaller dafa 
sef 

We demonsfrafe typical resulfs based on a choice of four radial bins. Each bin confains 90 
GRBs, wifhfhebins defined by 2.68 <z < 9.4 (group 1), 1.61 <z < 2.68 (group 2), 0.85 <z < 1.61 
(group 3), and 0 < z < 0.85 (group 4). Table 1 shows fhe significance of fhe null hypofhesis for 
fhis example using fhe Kolmogorov-Smirnov fesf fhaf fhe fwo disfribufions are differenf. Boldface 
type indicates fhaf fhe significance fhaf fhe fwo group’s 31sf nearesf-neighbor disfribufions differ by 
more fhan 3a. There are no significanf differences wifhin fhe group 1, group 3, or group 4 disfribu¬ 
fions, buf fhe 31sf nearesf-neighbor disfribufions in group 2 indicate a significanf anisofropy. The 
same indicafion of a large-scale anisofropy in fhis radial group is found from all nearesf neighbor 
disfribufions spanning fhe range 22 <k <55, indicafing fhaf fhe anisofropy occurs on an angular 
scale of infermediafe size. 

The significance, angular size, and locafion of fhe large, loose GRB cluster in fhe redshifl range 
1.6 < z < 2.1 can also be esfimafed using fhe boofsfrap poinf-radius mefhod described in section 
5 of Horvafh ef al. 2014 (0). This fesf compares fhe 44 GRBs found in fhe 1.6 < z < 2.1 radial 
bin fo 44 randomly selected GRBs drawn from fhe resf of fhe sample. This fechnique counfs fhe 
number of 1.6 < z < 2.1 bursfs wifhin a circle of predefined radius Q surrounding a random ’clusfer 
confer’ locafion. Comparison samples are created by randomly drawing 44 bursfs from fhe resf of 
fhe dafasef and counting fhe number of bursfs lying wifhin fhe same angular circle. Sfafislics are 
generated by repeating fhis process 10000 limes and counting fhe largesl number of bursfs found 
wifhin fhe circle during fhese runs. Once resulfs have been obfained, 0 is increased and fhe process 
is repealed for eighty differenf 6 values in equal area steps ranging from 0 = 12.84° fo 0 = 180°. 

The Monte-Carlo boofsfrap poinf-radius mefhod verifies fhaf fhe anisofropy found in fhe 1.6 < 
z < 2.1 redshifl range is consisfenl wifh angular clusfering. Forfy nine of fhe 80 angular cluslering 
scales tested exhibit excess numbers of bursts within the defined 0 limils. For example, using fhe 
measuremenfs from fhe cluster center locations producing fhe largesl GRB clusfers, an angular 
circle having a radius of 0 = 22.3° (corresponding fo 3.75% of fhe sky) is found lo confain 13 of 
fhe 44 bursfs (30%), a circle wifh radius 0 = 34.4° (corresponding lo 8.75% of fhe sky) confains 
18 of fhe 44 bursfs (41%), and a circle wifh radius 0 = 51.3° (corresponding lo 18.75% of fhe 
sky) confains 25 of fhe 44 bursfs (57%). Only fwo of 17500 boofsfrap cases had 25 or more GRBs 
inside fhis latter circle indicating a slalislically significanf (p=0.0001143) devialion (fhe binomial 
probability for fhis being random is = 2 x 10^^). 

The 42% increase in sample size should have noticeably decreased fhe significance of fhe 
1.6 < z < 2.1 clusfer if random sampling was responsible for ifs exisfence. However, fhe clusfer 
has become more pronounced in fhe 4° < 0 < 90° angular radius range as more GRBs have been 
added lo fhe sample. 

In Ibis range of angular radii, 49 angular circles confain enough GRBs lo exceed fhe 2a signif¬ 
icance level (compared fo 28 found in our previous analysis 0). Addifionally, fhere are 16 angular 
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circles containing enough GRBs to exceed the 3a level (compared to only 2 in our previously pub¬ 
lished result). Therefore, the evidence has strengthened that these bursts are mapping out some 
large-scale universal structure. 

Both of the aforementioned tests are independent of sky exposure; by calculating only the 
relative angular positions of the detected bursts the techniques assume simply that whatever biases 
present at different redshifts are the same as those in the 1.6 < z < 2.1 redshift range. 

3. Discussion 

Gamma-ray bursts are not distributed isotropically in the 1.6 < z < 2.1 redshift range; they ex¬ 
hibit evidence of large-scale clustering. This clustering, first identified in 2013 (||8|, ^]), has become 
more pronounced wifh recenf GRB defecfions by Swiff, supporting fhe idea fhaf fhe clusfering may 
be real rafher fhan due fo a sfafislical variafion in fhe defecfion rale. The nearesl neighbor lesf 
indicates lhal GRBs in Ibis redshifl range are likely fo have more neighbors af moderate angular 
separations fhan fhose al ofher redshifls. The fwo-dimensional poinl-radius mefhod also finds evi¬ 
dence for a large-scale angular clusfering in fhis redshifl range; fhe angular diameter encompassed 
by fhis clusfering is likely many lens of degrees across. 

The aforementioned techniques can be used fo demonsfrale fhe exisfence of a large, nebulous 
GRB clusfer, but the nature of the tests used here prevent us from knowing exactly where the cluster 
is located, what its structure is like, or how big it might be. Selection biases due to instrumental 
sky exposure and visual extinction by dust complicate this interpretation by reducing the rate of 
detection in some parts of the sky relative to others. Based on the detected bursts, the cluster covers 
roughly one-eighth of the sky, and seems to encompass half of the constellations of Bootes, Draco, 
and Lyra, and all of the constellations of Hercules and Corona Borealis. The name of the structure 
has been popularized as the Hercules-Corona Borealis Great Wall, or Her-CrB GW. However, 
we note that sampling biases could cause the cluster to be offset by many degrees from where it 
currently appears to be. 

Because GRBs are the most luminous, energetic objects known, they are tracers of the presence 
of normal matter that can be detected at distances where the matter is otherwise too faint to be 
observed. Based on the analysis performed here, we estimate the size of the Her-CrB GW to be 
about 2000-3000 Mpc in diameter. Few limits on its radial thickness exist, other than the fact that 
it appears to be confined fo fhe 1.6 < z < 2.1 redshifl range. This large size makes the structure 
inconsistent with current inflationary Universal models, as it is larger than the roughly 100 Mpc 
limit thought to signify the “End of Greatness” at which large-scale structure ceases. 

However, the Her-CrB GW is not the first optical/infrared structure found to exceed the 100 
Mpc size limit. In the 1980s, Geller and Huchra (||^) mapped galaxies and galaxy clusters in a 
portion of the sky to z ~ 0.03 and found a 200 Mpc size structure that was later called the CfA2 
Great Wall. In 2005 an object twice this size named the Sloan Great Wall ([^]) was reported. In the 
ensuing years, several other large filamentary structures have also been identified. Roger Clowes 
and his learn have found several large clusters of luminous quasars; fhe largesl of fhese being fhe 
Huge Large Quasar Group (Huge-LQG; [||]) having a lengfh of more fhan 1400 Mpc. 

As large as if appears fo be, fhe Her-CrB GW does nol necessarily violafe fhe basic assumptions 
of the cosmological principle (the assumptions of a homogeneous and isotropic universe). Theo- 
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retical large-scale structure models indicate that some structures will exceed the End of Greatness 
on purely statistical grounds (l|T^), and this may be one such structure (albeit a very large one). 
Along these lines, this may not be a single structure, but overlapping smaller adjacent and/or line- 
of-sight structures; the small number of bursts cuiTently found in the cluster limits our ability to 
angularly resolve it. In other words, this may become a semantic issue at some point, since a cluster 
of smaller structures might still be a larger structure. 

This research was supported by the Hungarian OTKA grant NNl 11016 and by NASA EP- 
SCoR grant NNX13AD28A. Discussions with E.G. Balazs are also acknowledged. 
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